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ABSTRACT 

We have obtained new resolved images of the well-studied HR 4796A dust ring at 18 and 25 /jm with the 8-meter Gemini telescopes. 
These images confirm the previously observed spatial extent seen in mid-IR, near-IR, and optical images of the source. We detect 
brightness and temperature asymmetries such that dust on the NE side is both brighter and warmer than dust in the SW. We show 
that models of so-called pericenter glow account for these asymmetries, thus both confirming and extending our previous analyses. 
In this scenario, the center of the dust ring is offset from the star due to gravitational perturbations of a body with an eccentric orbit 
that has induced a forced eccentricity on the dust particle orbits. Models with 2-yum silicate dust particles and a forced eccentricity of 
0.06 simultaneously fit the observations at both wavelengths. We also show that parameters used to characterize the thermal-emission 
properties of the disk can also account for the disk asymmetry observed in shorter-wavelength scattered-light images. 

Key words, circumstellar matter, planetary systems, Stars: individual: HR 4796A, Infrared: planetary systems 



1. Introduction 

HR 4796A is an A0V star with the highest fractional IR lumi- 
nosity (L/r/L*= 5 x 10~ 3 ) yet discovered a mong debris disks. Its 
recently revised distance estimate is 73 pc (Ivan Leeuweri 120071) 
(updated from 67 pc [Perryman et al. 1997]), a nd it has an M-sta r 
T Tauri companion with a separation of 7 .7" (lJura et al. [ [1991 . 
Its age of 8+2 Myr (IStauffer et all 1 19951) (based on the lithium 
abundance of HR 4796B) places the disk of HR 4796A in the 
very interesting early phase of chaotic disk evolution dominated 
by collisions (e.g., Kenyon & Bromley 2005). 

Based on the dust temperature, iJura et al ] <fl993h first tenta- 
tively suggested 40 AU (now 44 AU) as the main location of dust 
responsible for the HR 4796A excess thermal emission. Indeed, 
the mid-IR discovery images of the source at CTIO and Keck 
revealed a highly inclined ringlike d isk with a dust distr i bution 
peaki ng near 70 AU (now 76 AU) dJavawardhana et all [1991 
iKoerner et all Il998l) . Follow-up coronagraphic near-IR images 
with HS T/NICMOS tightly constrained the width of the ring to 
~17 AU (now ~18.5 AU), with the dust populati on severely de- 
pleted both interior and exterior to this region (ISchneider et all 
119991) . 

Further imaging observations of the thermal emission con- 
firme d a strong peak in the d ust density at 70 AU dTelesco et all 
2000; Wahhaj etafl 120051) . iTelesco et all (|2000) also noted a 
1.8-cr brightness asymmetry at 18 /im in which the northeast 
side of the disk is brighter tha n the southwest side. In a compan- 
ion paper, I Wvatt et all {1999) demonstrated that such an asym- 



metry could arise from the phenomenon of pericenter glow. A 
second body in the system on an eccentric orbit about the star 
may cause pericenter glow when, through secular perturbations, 
it effectively shifts the center of the dust ring away from the star 
and closer to the apastron of the perturbing companion's orbit. 
The side of the ring shifted cl oser to the star becom es relatively 
warmer and more luminous (iDermott et all 1998). giving rise 
to the "glow". At present, while several works have acknowl- 
edged the existence of a brightness asymmetry (e.g., Wahhaj et 
al. 2005, Debes et al. 2008), no alternative explanations for the 
brightness asymmetry have been proposed. 



We have obtained mid-IR images of HR 4796A with the 
Michelle and T-ReCS cameras at the Gemini North and Gemini 
South telescopes, respectively. These permit us to confirm un- 
ambiguously the brightness asymmetry, which was previously 
characterized at only the 1.8-cr level of significance. By imaging 
thermally emitting dust at two wavelengths, we also examine 
the spatial distribution of dust temperatures within the resolved 
disk, and we conclude that the dust population in the NE ansa 
of the disk is not only brighter, but also warmer, than the dust 
in the SW ansa. With these new data and with the consideration 
of shorter-wavelength images of scattered light, we re-examine 
the possibility that these observed asymmetric characteristics are 
caused by pericenter glow. 
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Fig. 1. Top: 18.1 -fim image contours drawn at 2-cr intervals start- 
ing at 5cr (background cr = 0.295 mJy/pixel). Bottom: 24.5-/im 
image contours drawn at l-cr intervals starting at 3-cr (back- 
ground cr = 1.463 mJy/pixel). Note that the 18.1-/im image was 
taken with Michelle (plate scale: 0.1 005 "/pixel) and the 24.5- 
fim image was taken with T-ReCS (plate scale: 0.089"/pixel). 
The images have been scaled to show the same area in square 
arcseconds, and both images are photosphere-subtracted. 

2. Mid-IR images of HR 4796A 

We observed HR 4796A with T-ReCS at Gemini South (using 
the Qb filter) in May 2004 and with Michelle at Gemini North 
(using the Qa filter) in April 2005. The central wavelengths for 
these two filters are 18.1 fim for Qa and 24.5 fim for Qb. A PSF 
and flux standard star were observed close in time to the target 
object for both sets of images. We reduced the images with an 
IDL reduction package developed at the University of Florida, 
and we performed photometric measurements with the Starlink 
GAIA program. The images in both bandpasses show that the 
disk of HR 4796A is clearly resolved, in the form of an elon- 
gated disk ~3" in extent with two prominent peaks in brightness 
(Figure [TJ. 

Table Q] lists results of the sky-subtracted aperture photom- 
etry of HR 4796A at 18.1 fim and 24.5 fan. The uncertain- 
ties reflect only measurement errors. Photometric uncertainties 
are driven primarily by variable sky transmission throughout 
the night, and lacking multiple standard star measurements, we 
adopt typical photometric uncertainties of 15% for both filters. 



Table 1. Photometry of the HR 4796 A Disk 



x c 


Total Flux Density 


Photosphere 


Disk Flux Density 


[jum] 


[mJy] 


[mJy] 


[mJy] 


18.1 


1106 ±7 


48 


1058 ± 7 


24.5 


3307 ± 47 


33 


3274 + 47 



Notes- Uncertainties listed in this table are measurement uncertainties. 
Photometric uncertainties are assumed to be 15% for both bandpasses. 



The photospheric contribution in the bandpasses of the imag- 
ing observation s are based on a 10.8-/vm value derived by 
lJura et all (1 1998b by extrapolating the g-ba nd (2.2 f im) flu x den- 
sity using a v 1 88 power law as given by the lKuruca (1 19791) model 
atmosphere with T = 9500 K and log g = 4.0. Photospheric val- 
ues at longer wavelengths were estimated with the assumption 



that HR 4796 A and Vega have the same 10-to-18 fim flux ratio, 
since they are both A0V stars. The photospheric contribution to 
the total flux density of the source is 4% at 18.3 fim and 1% at 
24.5 fim. 

The PSF was scaled such that the integrated flux matched 
the photosphere level, and its center was shifted to the posi- 
tion assumed for the star in the target source images. The max- 
imum intensities of the scaled PSF and of the disk (prior to the 
photosphere subtraction) at the stellar position are, respectively, 
1.4 mJy/pixel and 5.0 mJy/pixel at 18.1 fim and 0.3 mJy/pixel 
and 9.8 mJy/pixel at 24.5 fim. For the 18.3-/vm image, the 
adopted stellar position was the centroid peak of the 2-pixel- 
smoothed central lobe of emission (which we assume arises from 
the presence of the star) between the NE and SW ansae, where 
the peaks were also determined from a centroid measurement 
of the 2-pixel-smoothed image. The stellar centroid peak posi- 
tion in the 18.1 -/mi image differs from the actual midpoint be- 
tween the NE and SW peaks of emission by <1 pixel, which is 
insignificant as this is approximately the level of certainty in the 
determination of the peak location. One pixel in the 18.1 -fim im- 
age corresponds to 0.1", or ~1 AU. For the 24.5-ji/m image, the 
stellar position is not obvious, and so the midpoint between the 
two centroid peaks of the NE and SW ansae was adopted as the 
stellar position. The distance from the NE peak to the SW peak 
in the 2-pixel-smoothed image is 17.3 pixels (1.74") at 18.1 fim 
and 18.2 pixels (1.62") at 24.5 fim. 

The scaled PSF image in each bandpass was then sub- 
tracted from the corresponding image of HR 4796A. The total 
photosphere-subtracted flux density of the source was measured 
at each wavelength, and the results are given in TableQ] 



3. Disk asymmetry 

The global dust color temperature is 102 + 11 K, based on 
the 18. l-jt/m/24.5/im flux density ratio determined from aperture 
photometry of the entire photosphere-subtracted disk. This tem- 
pe rature is consisten t with the initial color temperature estimated 
bv lJura et al J (1 1 9931) based on IRAS photometry (T~l 10 K). The 
color temperature of the dust as a function of disk radius was also 
calculated using the approach described below. 

The position angle (PA), east of north, is defined as the angle 
between the vertical image axis and the line connecting the two 
central peaks of emission of the two ansae: 28.1° in the 18.1-yi/m 
image and 29.9° in the 24.5-fim image. We estimate an uncer- 
tainty for the PA of 0.9 degrees, half the difference between the 
values measured in the stacked image at each wavelength. We 
adopt a PA value of 29°, the average of the two values. For con- 
venient display, the images of HR 4796A were rotated counter- 
clockwise by an angle of 61°, to orient the disk plane parallel to 
the x-axis. 

Quantifying the asymmetric structure- The asymmetric 
structure is quantified in both the observations and in the model 
distributions (® by finding the brightness peak for each ansal 
"lobe" and its offset from the star. We calculated the following 
values for each image: 

- D mean , the mean offset of the lobes' peak brightness from the 
center, in pixels 

- 7p^-, the difference between the distance of the lobe bright- 
ness peaks (Dne - AswX divided by the mean offset of the 
brightness peaks from the center, given as a percentage value 

- F mea „, the mean peak brightness of the lobes 
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Fig. 2. Color temperature profile of the HR 4796A disk based 
on photometric measurements at 18.1 /mi and 24.5 /mi (images 
"cross-convolved" to achieve the same spatial resolution). The 
NE profile is represented by circles, and the SW profile is rep- 
resented by squares. Combined measurement uncertainties from 
both bandpasses are smaller than the size of the data points. The 
temperatures at the location of the peak dust density in the ring 
(labeled) are the key values to note. 




Fig. 3. Brightness profiles of the HR 4796A observations and 
model at 18.1 fim (lower curves) and 24.5 /mi (upper curves). 
The brightness profile for each observation or disk model was 
generated by summing the central 3" along the long axis of the 
disk. Model curves are shown for four different angles of peri- 
center, and we have concluded that a pericenter angle of zero 
provides the best fit. 



- yZ— , the difference in lobe peak brightness, (F^e - Fs w), di- 
vided by the mean peak brightness of the lobe (F mea „), given 
as a percentage value 

The level of the brightness asymmetry, AF/F mea „, is of pri- 
mary interest, and it is 15.3 + 2.6% at 18.1 /mi and 13.0 + 3.8% 
at 24.5 //m. The uncertainty in AF/F mea „ was calculated by tak- 
ing a circular symmetrical fit to the observations (i.e., repeating 
the modeling but imposing zero forced eccentricity) and then 
using Monte Carlo methods to add noise to the model at the 
observed level (0.007 mJy/pixel at 18.3 /mi, 0.135 mJy/pixel at 
24.5 /mi) to work out the asymmetry that would have been de- 
tected purely due to noise. 



Assessing the temperature asymmetry- After removal of the 
photospheric contribution, the angular resolution of each image 
was degraded to achieve the same resolution in both images for 
accurate spatial comparison in the color temperature calcula- 
tion. Using the gauss routine in IRAF, each of the target im- 
ages was convolved with a Gaussian profile having the same 
FWHM as the PSF from the other bandpass. The PSF FWHM 
was estimated by a Gaussian fit to the azimuthally averaged pro- 
file. The FWHM values were 5.19 pixels (0.52") in the 18.1-/mi 
image and 8.08 pixels (0.72") in the 24.5-/mi image. The PSF 
images themselves were convolved in the same manner, which 
confirmed that the resulting resolution (assessed with FWHM 
measurements) was identical for both images. The final resolu- 
tion of the cross-convolved images is ~0.9", which corresponds 
to ~65 AU for the source distance of 73 pc. 

The rotated images were cropped to a swath of vertical width 
Ay~3". The pixel values in these images were summed along 
the y-axis, resulting in a one-dimensional brightness profile for 
the disk at each wavelength. From these brightness profiles, a 
color profile and a color temperature profile (Figure|2|i were con- 
structed to illustrate variations in temperature along the extent of 
the disk. As a test of the robustness of the temperature estimates, 
the 24.5-/mi brightness profile was shifted by +1 pixel (~6 AU) 
along the x-axis relative to the 18.1 -/mi brightness profile be- 
fore the two brightness profiles were combined to determine the 
color profile. The two resulting color profiles that incorporated 
these offsets yielded color temperature profiles that did not differ 
from the results presented in Figure [2] within 60 AU of the star 
by more than 2 K. 

The NE and SW profiles are overplotted on the same radius 
scale (along the x-axis) to better show differences in tempera- 
ture at the same distances from the central star. We see no sig- 
nificant variation in color temperature within a radius of 15 AU 
from the star. While that result is expected for such a narrow 
annulus, this result is in fact attributable principally to the degra- 
dation of resolution that we carried out to make the tempera- 
ture estimates. This effect is also seen in the separation of the 
two ansae peaks, which decreased for both images following the 
cross-convolution, from 1.74" to 1.1" at 18.1 /mi, and from 1.62" 
to 1.58" at 24.5 /mi (where the lesser effect for the 24.5 //m im- 
age is likely due to being convolved by a smaller kernel). 

There is a clear temperature asymmetry, with the dust in the 
NE side of the disk having significantly higher color tempera- 
tures at greater than 30 AU from the star and in particular at the 
location of the emission peaks corresponding to the two ansae 
of the dust ring. The temperature difference peaks at ~9 K near 
60 AU. The uncertainties due to measurement error in the tem- 
perature estimates are less than 2 K. The color temperature un- 
certainties that incorporate photometric uncertainties are +8 K, 
but in this analysis we are most interested in the relative temper- 
atures of the two ansae and are therefore most concerned with 
uncertainties within the images and not absolute photometric un- 
certainties. 

4. Pericenter glow modeling 

4. 1 . Pericenter glow scenario and prior work 

One motivation for observing debris disks is the potential to re- 
veal as-yet unseen planetary companions through the disk mor- 
phology. A massive orbiting body can gravitationally perturb 
the orbits of smaller bodies like dust particles over long-period 
("secular") timescales, and these perturbations may be observed 
in the global characteristics of the disk. For example, a planet 
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on an eccentric orbit eventually imposes a forced eccentricity on 
the dust particles, which results in an offset of the disk's center 
of symmetry from the host star. While this offset may be ob- 
servable directly, other manifestations of the offset may be more 
prominent. In particular, if the dust population is azimuthally 
homogenous in size distribution and composition, the side of the 
disk (the pericenter side) that is offset toward the star will ex- 
perience enhanced heating, increasing both its temperature and 
brightness. The resulting asymmetry in brightness between the 
apocenter an d pericenter sides of the disk is referred to as peri- 
center glow (Dermo tt et al.L 1 1998b . 

The asymmetric morphology of the HR 4796A disk observed 
in previous thermal emission images at 18 pm was well ap- 
proximated by models th at invoked pericenter glow (IWvatt et all 
ll999UTelesco et alll2000h . but the statistical significance of the 
asymmetry in those earlier images was relatively low. We note 
that the b rightness asymmetr y was also previously detected at 
10.8 pm (Telesco et al., 2000), but the much greater prominence 
of the starlight at those wavelengths makes such shorter wave- 
length observations less useful for our analysis. With new 18.1- 
yum images and the addition of images at 24.5 /mi in this dataset, 
we have re-examined this hypothesis by using the same models 
to reconstruct the observed thermal emission and color tempera- 
ture profiles. 

4.2. Thermal emission reconstruction 

Model inputs- The physics of pericente r glow and the assoc iated 
modeling of it are outlined in detail bv lWvatt et al.l (11999). and 
here we review only the key parameters involved in that analysis. 
The model for the disk's density distribution of dust was gener- 
ated with three distinct inputs: the physical structure of the disk, 
the combination of the optical properties and size distribution 
of the dust particle population, and the orientation of the disk 
along the line of sight. We recall that the mass of the perturbing 
body itself does not come into play and is not constrained by 
this approach. The physical structure of the disk is modelled by 
a radial density distribution of dust with defining parameters as 
follows: a m i„ and a max , the inner and outer semimajor axes of the 
dust ring (which are roughly equivalent to the inner and outer 
radii), y, the power-law index of the semimajor axis distribution, 
<?/, the forced eccentricity, and cr, ol , the total surface area of the 
dust particles in the disk. The forced eccentricity is that imposed 
on all particles in the disk by the planet's secular perturbations 
dMurrav & Dermotlll999l) : as such, the planet's orbit is not af- 
fected by its own perturbations. The eccentricities of particle or- 
bits are the vector sum of this forced eccentricity and their own 
proper eccentricity. To first order, the disk remains circular but is 
just offset. The total surface area of the dust is effectively a scal- 
ing factor that can be adjusted to approximate the overall bright- 
ness level of the dust emission. While the intrinsic disk offset 
and asymmetry (i.e., if viewed face-on) are determined by the 
forced eccentricity e/, the measured offset and asymmetry are 
determined by the additional variables of disk inclination, i, and 
the angle of pericenter, (Of. For example, the forced eccentricity 
required to replicate a given measured asymmetry increases as 
the angle of pericenter deviates from being perpendicular to the 
line of sight. 

The optical properties of dust were calculated us- 
ing Mie theory, Rayle ig h-Gan s theory, and g eo metric 
optics (|Li & Greenbergl 1 1997b lAugereau et all Il999t 
Bohren & Huffman, Il983l) for astronomical silicate spheres the 
diameter of which (D np -2 pm) was constrained by a simul- 
taneous fit to the images at both wavelengths. The resulting 



dielectric constant was 6.7xl0~ 5 , and the absorption coefficients 
were 0.54 and 0.30 at 18.1 and 24.5 pm, respectively. While a 
single particle size is used to represent a particle population that 
has a range of sizes and composition, it reasonably encapsulates 
key mid-IR emission characteristics of the disk. We adopted 
values for the stellar luminosity and effective temperature of 
21 L and 9500 K, respectively (e.g., Debes et al., 2008). 

The output from the disk radial distribution model and the 
calculation of the dust optical properties were combined with 
the specified orientation to calculate the line-of-sight brightness 
and generate images of the model disk. The model images at 
each wavelength were rotated to the same position angle as the 
observed disk, and each image was convolved by a PSF kernel 
with the profile width of the observed PSF reference star in the 
same filter. In the generation of color temperature profiles for the 
model data, the images were also convolved by the profile width 
corresponding to the other filter to match the treatment of the 
observational data ($3). 




Fig. 4. Color temperature profiles of the HR 4796A observations 
and models. Colors were calculated by summing the flux within 
a 3"-wide swath along the long axis of the disk for each image. 
Model curves are shown for two different angles of pericenter, 
and we have concluded that a pericenter angle of zero provides 
the best fit. 



Optimizing the model- With the images of the model disk, 
we generated brightness profiles and color temperature profiles 
in the same way that we did with the observational data (§01. 
We compared the observed and model profiles, and we adjusted 
the model parameters until the model disk yielded a good ap- 
proximation to the observations. First, we varied the inner and 
outer disk radii and the total disk surface area until the brightness 
profile width and magnitude were well approximated. We then 
adjusted the forced eccentricity until the brightness asymmetry 
level was matched for a fixed value of the angle of pericenter. 
The fit was repeated for several values of the angle of pericen- 
ter to assess the corresponding forced eccentricity necessary to 
reproduce the same asymmetry (Table |2J. The forced eccentric- 
ity and the angle of pericenter define the magnitude of the ra- 
dial offset asymmetry in the disk model, which in turn defines 
the brightness asymmetry and the dust color temperature asym- 
metry. The initial assessment of how changing these parameters 
affected the resulting brightness profiles was performed by chi- 
squared minimization on the brightness profiles and images. The 
brightness profiles of the observations and the model disk with 
varying values for the angle of pericenter are shown in Figure[3] 
and the color temperature profiles are shown in Figure [4] The 
quantitative measurements used to assess the asymmetry in the 
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Table 2. Asymmetry parameters for observations and models 



18.1 24.5 







e f 


u mean 


[%] 

t mem 


■£t [%] 


3^- [*] 


reduced 


observed 









15.28 





13.00 




model a>f 


= 0" 


0.06 


-11.71 


15.26 


-6.11 


13.10 


2.17 


model a>f 


= 40° 


0.10 


-10.81 


15.31 


-7.01 


13.30 


2.31 


model ojf 


= 55° 


0.13 


-11.18 


15.35 


-7.15 


13.05 


2.21 


model ojf 


= 75° 


0.30 


-14.50 


10.55 


-7.39 


12.70 


2.26 



Notes. All models were simultaneously fitted to 18.1- and 24.5-//m brightness profile linecuts of 0.5" and share the following parameter values: 
a m i„=lQ AU, £Wr = 84 AU, y=-1.5, D lyp =2 yum, and i=14.1°. Also recall that the value j^— for the observations is because we defined the stellar 
position at 24.5 yum as the midpoint between the brightness peaks of the two ansae (with the measured difference between the midpoint and the 
stellar peak at 18.1 prntobe <1 pixel.) 



model and in the observed disk, as described above, are summa- 
rized in Table [2] 

Results of the modelling- We have determined that a forced 
eccentricity of 0.06 imposed on a ring of dust with D typ - 2 ftm 
spanning 74-84 AU can reproduce the thermal IR images, fit- 
ting the observed brightness asymmetry well at both 18.1 and 
24.5 yum. The disk density power-law exponent was set to -1.5. 
The same result for the characteristic particle size was inferred 
previously in order to simultaneously fit 10-/um and 1 8-//m im- 
ages of HR 4796A (IWvatt et all 1 19991: iTelesco et alll2000h . We 
have estimated the uncertainty in the forced eccentricity to be 
0.01, which is the change required to effect a l-cr change in 
the brightness asymmetry of the model at 18.3 fim (the more 
limiting case). The angle of pericenter was 0+30° (perpendic- 
ular to the line of sig ht), and the disk inclination was 14.1° 
dSchneider et all 2009h. We note th at the se results diff e r from 
those found by IWvatt et alj d 19991) and ITelesco et al] {2000) 
(e/=0.02, a)/ =75°) likely because of the different levels of ob- 
served asymmetry, which are nonetheless within 3cr of one an- 
other (5.1 + 3.2% in the previous work, 15.3 + 2.6% in this 
work). 

We have also used the model to quantify the limits set by 
these data on the presence of hot dust inside the ring. We com- 
pared the photometry of the photosphere-subtracted images with 
that of the model within an aperture radius equal to the size of the 
PSF FWHM (0.52" at 18.3 //m, 0.72" at 24.5 yum), and then we 
repeated the measurement with the same aperture on the models, 
where an unresolved flux component could be added at the po- 
sition of the star. Excluding calibration uncertainties, the mea- 
sured flux density within a 0.52"-radius aperture on the 18.3- 
yum photosphere-subtracted image was 289 + 3 mJy. The same 
measurement on the model with no added unresolved flux com- 
ponent was 290 mJy. Therefore, the maximum unresolved con- 
tribution we can add to the model that is 3-cr consistent with 
the observations is 8 mJy. Repeating the same procedure at 24.5 
um with an aperture of radius 0.72", we measured a flux den- 
sity of 1290 + 12 mJy for the photosphere-subtracted images 
and 1289 mJy for the model. Therefore, the 3-cr upper limit of 
1326 mJy for this aperture allows the presence of an unresolved 
component at a level of 37 mJy. However, these measurements 
do not account for the assumed 15% calibration uncertainty. 
Since a change in calibration would also result in a correspond- 
ing change to the model (to fit the peaks at the correct level), it is 
instructive to consider what would have been derived with a re- 
vised calibration factor. In the extreme situation that the calibra- 
tion factor was 45% higher, the flux measured in the previously 
described apertures at 18.1 and 24.5 yum respectively would have 
been 438 + 4 mJy and 1881 ± 1 7 mJy, with corresponding model 
fluxes of 290*1.45 = 421 mJy and 1289*1.45 = 1869 mJy, and 



so a flux excess of 17 + 4 and 12+17 mJy. Therefore, a 3-cr 
deviation in both calibration and statistical uncertainty puts the 
limits on the u nresolved flux at < 29 mJy and <63 mJy at 18.1 
and 24.5 yum. Wahhaj_etalJ([2005) estimated an exozodiacal con- 
tribution of 87 ± 22 mJy at 24.5 yum, which is consistent to within 
3cr with the upper limit we derive. 

4.3. Scattered light comparison 

A brightness asymmetry has also been detected at wavelengths 
that are do minated by light scatte red by the dust disk. Images 
obtained bv lSchneider et alJ (Tl999) at 1.1 /im were suggestive of 
a brightness imbalance favoring the NE that w as later noted to 
be at the 10-15% level dSchneider et all 1200 ll) . The disk ansae 
were again compared in the analysis of data sets taken with HS T 
NICMOS at 1.71-2.22 fjxa, but the NE ansa was n ot consis- 
tently brighter in all bandpasses dDebes et all 120081) . Most re- 
cently, HST STIS images at 0.57 fim show that the SW ansa 
is 0.74 + 0.07 times as bright as the NE ansa, as measured in 
"lobes" within +20° of the major axis and between 0.8" and 1.5" 
from the star. ISchneider et a"LT(l2009l) measured an offset asym- 
metry, {Dne - Dsw)/D mecm = -3.7%, such that the center of the 
disk is offset from the star by 1.4 + 0.4 AU (19 + 6 mas) in the 
plane perpendicular to the line of sight. 

Thus, another test for our pericenter glow model is whether it 
can replicate the observations in scattered light. To examine this, 
we modified the model's particle properties to correctly calcu- 
late the scattered light fluxes, using a combination of Mie theory, 
the Rayleigh-Gans approximation, and geometric optics, which 
yielded a mean albedo of 0.05. We assumed that the same parti- 
cle population (in this case, modeled as a first approximation by 
2-//m silicate spheres) is responsible for both the thermal emis- 
sion and the scattered light reflection. We adopted a Henyey- 
Greenstein phase function to model the non-isotropic scattering 
properties of the small grains (lAugereau et all 1 1 9991) . with the 
form 



/(«) 



1 



1 - g 2 - 2gcos(a) 



(1) 



where a is the scattering angle to the line of sight and g is 
the asymmetry parameter. The parameter g indicates the rela- 
tive level of front or back scattering, such that g = 1 for 100% 
forward scattering, g = — 1 for 100% backward scattering, and 
g = for isotropic scattering. We adopted th e value g - 0.16 as 
determined from the o ptical S TIS images dDebes et all [2008; 
ISchneider et all 120091) . wh ich is consistent w ithin 2<x of the 
values for g determined by iDebes et"al] ((2008) for the near-IR 
NICMOS images. 
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To check that our model produce d realistic flux levels, we 
compared them to those measured by Schn eider et al.1 (120091) at 
a wavelength of 0.57 pm. Due to the presence of coronagraphic 
spikes in the image, the flux density measurement had to be cor- 
rected for the missing area, yielding a total flux 9.4±0.8 mJy. The 
total flux in our model scattered light image is 10.8 mJy. As a 
secon dary check, we used the formula given bv lWeinberger et all 
(1999) to predict the surface brightness of a disk in scattered 
light (in mJy arcsec~ 2 ) 



where F is the flux received from the star in mJy, r is the optical 
depth of the scattering material, a> is the albedo (from Mie the- 
ory, where {2 sra =0.08), and the factor tcj can be approximated 
by Ldisk/Li, if co — 1. The expected flux density calculated with 
this relation agrees with the flux density values both from ob- 
servations and from our model to within 20% at 0.57, 1.1, and 
1.6 yum. 

We also replicated the lobe asymmetry measurement by 
ISchneider et alj d2009l) with our disk model. Within the same 
apertures (within ±20° of the major axis and between 0.8" and 
1.5" from the star), we determine a SW/NE brightness ratio of 
0.92 at 0.57 fim and 0.93 at 1.1 /m i. Our model values are there- 
fore consistent within 3<x with the Schneid er et alJ ([2009) value 
of 0.74 + 0.07, and are also consistent with a simple 1/r 2 at- 
tenuation of starlight due to the disk ansae being at different 
distances from the star. Our model also replicates the radial off- 
set asymmetry at a lev el of 3.3%, which is consistent with the 
ISchneider etail (120091) val ue of 3.7 + 1.1%. 

ISchneider et al.l (120091) acknowledge that dynamical pertur- 
bations such as those we propose may be responsible for part 
of the brightness asymmetry, also suggesting that a segregation 
of the particle population or an azimuthal density variation in a 
homogenous dust ring could produce the asymmetry. However, 
whether such a model could also reproduce the observed tem- 
perature asymmetry has yet to be tested. Additionally, we note 
that particles in the disk should be moving slightly faster at 
pericenter, relative to their velocity at apocenter, and the colli- 
sional timescale at pericenter should be relatively shorter in turn. 
Therefore, the population of small particles that scatter near- 
IR light may be enhanced due to the increased collision rate 
near pericenter, thus yielding a density variation in the dust ring 
and a greater brightness imbalance in scattered light. There is 
the potential to further develop the pericenter glow model with 
such an enhanced particle population, which may help to deter- 
mine whether factors beyond dynamical perturbations must be 
invoked to explain the observed morphology. Such refinements 
to our model would bear greater consideration especially if the 
large near-IR brightness asymmetry is confirmed. 



assist in confirming whether the pericenter glow model is suf- 
ficient to explain the disk asymmetry observed throughout the 
spectrum. 

The model disk is comprised of particles 2 fim in diam- 
eter (D tVD ), i n agr e ement with previou s analysis of the disk 
dWvatt et all 119991; iTelesco et all l2000h , and the particles are 
assumed to be spherical silicates. Dust particles smaller than 
~8 //m in diameter should be blown out roughly within an or- 
bital timescale (~400 yr). In reality the dust population is not this 
simple, but we have adopted this value to characterize the overall 
population. We must consider that the particles may have a more 
complex set of optical properties that allows them to be heated 
to the temperatures we infer and to remain in the system without 
being rapidly blown out. More complex dust properties (e.g., Li 
& Lunine 2003) could be incorporated into the pericenter glow 
model in future iterations. 

Finally, we recall that a key result for a system like 
HR 4796A that exhibits pericenter glow is the implied presence 
of a perturbing companion. In fact, the principal perturbing body 
may be the nearby M-star companion HR 4796B. However, the 
orbital parameters of HR 4796B are unknown and may not be 
commensurate with its being the primary perturber. The sharp 
inner edge of the ring points to the possibility that the perturb- 
ing object responsible for the asymmetry is also tru ncating the 
inner edge, as has b een suggested for Fomalhaut (iKalas et all 
2005: lOuillenl 120061) . In this case, the putative planet responsi- 
ble would be expected to lie just interior to the ring. For example, 
a Jupiter-mass planet would be expected to orbit at ~60 AU with 
an eccentricity of ~0.06 (maximum projected separation 0.87"), 
but it could be closer to the dust ring if it were more massive. 
Given the ambiguity in potential perturbers and the fact that the 
level of forced eccentricity induced by the companion does not 
constrain its mass, the continued search for the putative compan- 
ion is warranted, if not with the current generation of telescopes 
then with the next, such as the s pace-based JWS T or the ground- 
based TMT or E-ELT facilities. IKalas et alj (120081) have demon- 
strated a successful direct planet detection that was pursued fol- 
lowing similar evidence of a forced stellar offset, as in the case 
of Fomalhaut, which is nearer and brighter than HR 4796A. 
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5. Discussion and conclusions 

Using new mid-IR images at 1 8 and 25 fim, we have confirmed 
the previously obse rved brightness asymmetry discovered by 
ITelesco et ail (l2000l) in the HR 4796A disk to a high level of 
statistical significance. We have determined that a model disk 
with forced eccentricity of 0.06 that is dynamically induced by a 
companion can replicate that brightness asymmetry. After incor- 
porating the dust scattering properties, we can also reproduce, 
to within 3<r, the level of brightness asymmetry seen in near-IR 
scattered-light images . We believe that near-IR imaging in ad- 
ditional bandpasses will further constrain the level of the asym- 
metry observed in scattered light, and that such constraints will 
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